Critical Solutions of Three Vortex Motion in the Parabolic Case 
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^ ' Abstract 

Grobli (1877) laid the foundation for the analysis of the motion of three point vortices in 
a plane by deriving governing equations for triangular configuration of the vortices. Synge 
(1949) took this formulation one step further to that of a similar triangle of unit perimeter, 
via trilinear coordinates. The final reduced problem is governed by an integrable two- 
dimensional system of differential equations with solutions represented as planar trajectories. 
Another key to Synge's analysis was his classification of the problem into three distinct 
cases: elliptic, hyperbolic and parabolic corresponding, respectively, to the sum of products 
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of pairs of vortex strengths being positive, negative or zero. The reduction of the vortex 
configuration, a curve in space to a planar curve is one-to-one, except along a critical planar 
curve C in the parabolic case. Each point on C represents a triangle of unit perimeter 
corresponding to a family of similar vortex configurations, expanding or contracting. The 
latter would lead to coelescence of the three vortices. Tavantzis and Ting (1988) filled most of 
the gaps left by Synge regarding the dynamics of the problem, and showed in particular that 
points on C corresponding to similar expanding families of vortex configurations are stable 
while those corresponding to similar contracting families are unstable. Their investigations 
yielded an exhaustive description of the motion and stability of three vortices in a plane 
except for the global behavior of the vortex configurations in a narrow strip containing C. 
The main contribution of this paper is a complete description of the global dynamics in such 
a strip, which cmphaticalUy demonstrates that three distinct vortices almost never coalesce. 

Keywords: Three point vortices; Trilinear coordinates, Parabolic case; Critical curve. 



1 Introduction 

The planar motion of three point vortices in an incompressible fluid was studied by Grobli in 
1877 [Ij and Synge in 1949 [2] among others. With the aid of the integral invariants, Grobli 
uncoupled the autonomous six-degree of freedom problem to the study the vortex configura- 
tion and the translation and rotation of the configuration and essentially demonstrated the 
integrability of the system. The vortex configuration problem refers to the deformations of 
the triangle, Ap{t), formed by the three point vortices at Zj{t) in the complex z-plane. The 
triangle is defined by the lengths, Rj{t), j = 1,2,3, of its sides, with the j-th side facing the 
j-th vortex. The subscript p denotes the perimeter, p = Ri -^ R2 + R3, shown in Fig. [T] (a). 
It is a three degree of freedom problem and is integrable, with the integral curves, or spatial 



trajectories, R(t), in Cartesian coordinates, Rj{t), given by Synge [2]. The solution was further 
symphfied by Synge [2j who introduced the trihnear coordinates Xj, j = 1, 2, 3 with Xj = Rj/p- 
As shown in Fig. [1] (b), the spatial integral curve, R(t) is projected radially onto the plane, V, 
which intersects the j-th axis at Pj with intercept \/2/3. Its cross-section in the first octant 
is the equilateral AP1P2P3 with side 2/v3 and height 1, as shown in Fig. [1] (b). The heavy 
dot on the plane V with position vector ^2/3R/p, denotes the radial projection of R. The 
trilinear coordinates, rcj's, in turn denote the distances from the heavy dot to the sides of the 
triangle P1P2P3 opposite to the vertices -Pj's, as shown by the thin lines. Note that the sum, 
xi + X2 + X3 = 1, is the height of the triangle PiP2^3- The trilinear coordinates represent the 
sides of Ai(t) with perimeter 1, which is similar to and has the same orientation as Ap in Fig. [1] 
(a). 
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Figure 1: (a) The triangle Ap formed by the three point vortices at Zj{t), j = 1,2,3, in the 
complex z-plane, and (b) the coordinate axes, i?i,i?2 and R3, the plane V or the q/3 plane and 
the trilinear coordinates Xj, shown as the lines from a point, the heavy dot, on V to its edges 
in the first octant. 



Because of the triangle inequality, the admissible solutions, Xj of Ai, are confined to 
^QiQ2Q3 and a point on its positive (negative) side with normal vector pointing away from 



(towards) the origin O corresponds to a configuration Ap with the three vortices in counter- 
clockwise (clockwise) orientation. Two points on the opposite sides of V with the same Xj are 
called the image points of each other. 

Synge p] identified the sum of the products of the vortex strengths, kj, 

K = kik2 + k2hi + k^ki, (1) 

as the primary parameter, and described the topologies of the integral curves and critical points 
in the plane V according to the classifications, 

elliptic; i^ > 0; parabolic; K = Q] and hyperbolic, K < 0. (2) 

In particular, he studied the critical point, E, at the centroid of APiP2^3 and also that of 
l\QiQ2Q'i (Aq for short hereafter) on the positive side of 7-", where Xj = 1/3. The corre- 
sponding configuration is an equilateral triangle and was shown to be stable (unstable) for 
the elliptical (hyperbolic) cases. This is true also for the image of E, the centroid E* on the 
negative side of V. 

For the parabolic case, Synge pointed out that the critical point E lies on a critical curve 
C and partitions C into expanding and contracting branches, while the images of those two 
branches on C* are contracting and expanding, respectively. The study of the local stability 
and global trajectories near these two branches on C and their images is the main subject of 
this paper. 

The planar vortex problem was studied via Hamiltonian formalism by Lin in 1943, [lOj, 
Novikov in 1975 [3], Aref in 1979 @] among others. Tavantzis and Ting in 1988 [5] continued 
Synge's analysis via trilinear coordinates. For each classification by X in ([2D, they located all 
the critical points and studied their stability and described the topologies of the integral curves 
and separatrices in the plane V- 



Numerical studies of the motions of three vortices in a half-plane were carried out by Knio et 
al. ^J. for vortices of equal strength, k, i. e., K = 3A;^ > in the elliptic cases. The numerical 
results showed that the motion of the vortices can be regular or chaotic depending on their 
initial configuration, which was taken to be either an equilateral triangle or collinear. The 
transition from regular to chaotic for those half plane problems can be anticipated when the 
distances between the vortices are much smaller than their distances from the boundary. The 
half plane problems can then be considered as perturbations of three vortices in the entire plane, 
for which the equilateral configuration corresponds to a center while the collinear configuration 
corresponds to a saddle point in elliptic cases, |2]. This observation motivated a sequence 
of theoretical and numerical studies of perturbed three vortex problems in parallel with the 
classical perturbed two body problems. Here the perturbations arise for example from the 
presence of a boundary and/or additional vortices far away from the three main vortices or 
from the approximations the interactions of three coaxial vortex rings by three vortices in a 
meridian plane when the distances between the rings are much smaller than the radii of the rings. 
See the article by Blackmore et al. [7] and the references therein and recent studies by Ting et 
al. [8j and Knio et al. [9j. Several numerical examples were presented to show transitions from 
periodic and quasi-periodic to chaotic regimes in accordance with the theoretical results. 

This paper does not deal with perturbations. Rather, it presents an exhaustive study of 
the global behavior of the dynamics of three vortices in a plane in parabolic cases, with initial 
vortex configurations in a narrow strip containing the critical curve C. This completes the local 
stability analysis of C presented Tavantzis and Ting [5] . 

For the strengths of three vortices, at least two of them have to have the same sign. Thus 
we can always choose the sign of those two to be positive and denote them by ki and k2 with 

ki>k2> 0, and k^ = -kik2/{ki + /C2) < 0, (3) 



for the parabolic case, K = 0. After ordering the three vortices by Eq. ([3]), the initial positions of 
the vortices Zj{0) in the z-plane define the sides Rj{0) of the triangle Ap and hence the trilinear 
coordinates x(0) in the plane V and in addition the orientations of the vortex configuration, 
clockwise or counterclockwise, denoted by the index 7 = +1 and —1, respectively. To account for 
the opposite orientations, 7 = ±1, Synge made use of both sides, the positive and negative sides, 
of V, or AP1P2P3. The normal vector to the positive or the negative side, points away from or 
towards the origin O, respectively, (see Fig.[T|). He assigned the points Xj with 7 = +1 (7 = —1) 
to be on the positive (negative) side of AP1P2P3. A change of orientation can take place only 
when crossing over an edge of the triangle. 

It was pointed out by Synge [2j that the mapping from a point R on the spatial integral curve 
to its radial projection onto a point on the planar curve in P, is one to one if the velocity R(t) 
is not in the radial direction, i. e., R x R 7^ 0, which can occur only on a planar critical curve C 
in a parabolic case. The critical curve is partitioned into expanding and contracting branches. 
Corresponding to a point on the curve, C, the sides Xj's of the Ai remain stationary, while Ap(t), 
similar to Ai, is expanding, with p > on one branch and contracting p < on the adjacent 
branch. It was shown in [5j that the similar expanding solution is stable while the contracting 
one is unstable. Similar contracting solutions, albeit unstable, have received attention because 
they could in principle lead to simultaneous colliding or coalescence of three vortices (see for 
example Newton [11], pp. 84-88). Thus we need to reinforce the local instability analysis of 
[5] , for a point on a contracting branch of C by nonlocal analysis leading to the conclusion that 
coalescence of three vortices is unlikely. 

In Sec.[2]we present a brief review of the three point vortex problem, define the symbols, and 
summarize the relevant results of Synge [2j and Tavantzis and Ting [5j. To render this paper 
self-contained, we reproduce derivations and descriptions of results relevant to the current 



investigation, namely, the spatial trajectories of vortex configurations and their reduction to 
planar trajectories via trilinear coordinates. 

In Sec. O we study in detail the planar trajectories and the critical curve C in parabolic 
cases. We extend the local stability analysis of the critical curve in f5] by showing that a small 
deviation off C from a point on a contracting branch produces a trajectory that departs from 
the neighborhood of x and is finally attracted to an expanding branch of the critical curve. 
Depending on the location of the point x on C and the perturbation of x to a point above or 
below C, we identify three types of trajectories. An outline of this investigation was presented 
recently by Ting and Blackmore |12j . 

In Sec. m we present numerical examples to demonstrate how the the configuration Ai 
departs from the neighborhood of a contracting branch of the critical curve C along a trajectory 
to the expanding branch. We show the three types of trajectoties and the variations of the 
configurations depending on the locations of initial point x, as predicted in Sec. [3l 

2 Background 

Adopting the symbols of [2] and [5] , we relate the strength of a vortex kj to its circulation Tj 
by kj = Tj/2'K, j = 1, 2, 3, and denote the position of the j-th vortex in the complex plane by 
Zj{t) = Rzj + iQzj, and the length of the side of Ap facing its j-th vertex at Zj by Rj. 

Without loss of generality, we assign the strengths of the three vortices according to Eq. ([3|) . 
We set the initial perimeter as the length scale, p{0) = 1 and then the time scale is defined by 
p^(0)//c2- The equations of motion of three point vortices in a planar incompressible inviscid 
flow are: 

ijit) = -iJ2^^^ for i = l,2,3. (4) 

The summation in m ranges over 1,2, and 3. We obtain from Eq. (j4j) the well-known theory of 



KirchhofF or the conservation laws of vorticity and the moment of vorticity in the z-plane ( [13] , 
pp. 229-230), 

2_]/sjij = 0, and 2_] ^i^i'^i(^) ~ ^^ (^) 

j J 

The second equation in polar coordinates yields 

y^ kjPj = const, and ''^kjPjOj = K, (6) 

j 3 

where pj and 6j denote the magnitude and argument of Zj, respectively. 

In the next section, we will reproduce the analyses of Grobli (1877) and Synge (1949) which 
decoupled the six equations for the real and imaginary parts of Eq. ([4]), j = 1,2,3, to three 
equations for Rj, the vortex configuration and three for the motion of the configuration. The 
equations for the latter will come from the above two conservation equations in ([5]). 

2.1 Deformation of the configuration Ap 

The deformations of Rfs are governed by linear combinations of Eq. ([4]). For example, from a 
linear combination of the first two equations dH) , we get 

(zi - Z2)(ii - i2) = i?3i?3 + iRli^2i = i{ki + ^2) - iksRslR^^e'^^^^-^^'^ - R^^e'^^^^~^^^^]. (7) 

Here ipji denotes the argument of the side zi — Zj, and ipij = iIjji + vr denotes that of Zj — zi. The 
real part of the equation yields the rate of change of the side -R3 = \zi — Z2\. Likewise, Synge 
[2] obtained the rates of change of all three sides of Ap, 

Ri ^ R2 ^ R3 ^ 2-f\A\ 

kiRi[Rr^ — ^2) k2R2\R\ — ^3) k^R3[R2 ~ Ri) R1R2R3 

where \A\ = [s{s — Ri){s — R2){s — R^)]^'"^, with s = p/2, stands for the area of Ap. From the 
initial data for the Zj's, we have the data for the i?j's of the A/ and then Eqs. ([8|) define the 
spatial integral curve, R(t), in the first octant, where Rj > 0, j = 1, 2, 3. The integral curve 7^ 



is defined directly by the two integral invariants, or by the intersection of two surfaces, 

k^^RJ + k^^Rl + k^'^RJ = initial value a (9) 

^i/fei^i/fe2^i/fc3 ^ initial value 6. (10) 

They were obtained by KirchhofF for A^-vortices, with A^ = 3 here, (see e. g. [13[ page 230]). 
From R(t), we define the direction along the integral curve and call the directed integral curve 
the trajectory of the configuration, or simply the trajectory. 

To recover the positions Zj (t) corresponding to a point on the integral curve TZ, we need to 
find three combinations of the primary six equations in Q independent of the three equations 
for the configuration ([8]). These three combinations appear readily from the real and imaginary 
parts of the conservation of the center of vorticity in ([5]) and the imaginary part of the polar 
moment of vorticity in Eqs. ([6]). The real part of the equation for the polar moment is not 
independent because it is equivalent to Eq. ([9]), see [2j. 

We first find the dependence of Rj and their inclinations, say R3 and ^12, on t by integrat- 
ing Eq. ([8]). That would then define the temporal variations of the configuration Ap and its 
inclination. Finally the locations Zj{t) of the vertices of Ap(t) are defined by the first invari- 
ant in dS]), or by locating the (stationary) weighted center of Ap(t), Zc = (A;izi(0) -|- k2Z2{0) + 
k3Z3{0))/{ki + k2 + k^), for ^1 + ^2 + ^:3/ 0, which holds for parabolic cases. 

The reduction of the spatial curve R to a planar curve by Synge [^ via trilinear coordinates 
is outlined in the next subsection. 

2.2 The trilinear coordinates 

Synge [2] projected the point R radially onto the point, q on the plane V, 

Ri + R2 + R3 = a/s/S, (11) 



which intersects the axes respectively at Pi, P2 and P3, with equal intercept, y^2/3. The plane 
V intersects the coordinate planes at the same angle, arcsin \/2/3. As shown in Fig. [1] (b), 
the equilateral AP1P2P3 denotes the section of V in the first octant with each side 2/\/3 and 
height H = 1. The trilinear coordinates Xj, j = 1, 2, 3 of point q in P are the signed distances 
from the sides of AP1P2P3. The sign of the distance from a side is positive (negative) if the 
distance from the side to q points inwards to (outward from) AP1P2P3. The sum of the trilinear 
coordinates is always equal to the height H, 

Xi + 2;2+X3 = l. (12) 

The Xj's are all positive when q lies inside AP1P2-P3 and are related to the Cartesian coordinates 
R'j of q by, 

— = V -, while -f = — ^ and -^ = p, j = 1,2,3. (13) 

Xj V o JAjj "y A I o Xj 

Thus the Xj's also denote the sides of Ai with perimeter 1, similar to Ap with the same 
orientation. The mapping of the spatial trajectory R(t) to the planar trajectory x(i) is one-to- 
one, provided that the radial projection of R onto V is nonzero, or the spatial trajectory is not 
radial 

R X R / 0. (14) 

This condition is fulfilled for the elliptic and hyperbolic cases. The condition can be violated 
only for the parabolic case, where the spatial trajectory can be radial. Then the radial projection 
of the trajectory onto V reduces to one singular point q. The loci of these singular points is 
called the critical curve C such that each point q on C represents a Ai corresponding to a radial 
trajectory of Ap r^ Ai either moving away from the origin with p > or moving inward with 
p < 0, see |2j and |5j. More details on the critical curve will be elaborated later in Sec. [3l to 
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set the stage for the mam objective of the current paper, which is to identify branches of the 
critical curve that attract or repel nearby trajectories. 
Due to the triangle inequality, 

1 1 1 

Xl+X2>->Xz, X2 + X^>->Xl, X^ + Xl>->X2, (15) 

X has to lie in the Ag with vertices Qj, j = 1,2,3 lying at the midpoint of the j-th side of 
APiP2-P3- The equality signs in ([15]) hold on the edges Q1Q2 , Q2Q3 and Q3Q1 of Aq, where 
X3, xi and X2 = 1/2, respectively, and the three vortices are collinear with the area of Ap = 0. 
The vertex Qj of Ag corresponds to a singular case where the three vortices degenerate to two. 
For example, at Qi, where xi = 0, and X2 = X3 = 1/2, ov Ri = and R2 = R3, the three 
vortices degenerate to two, one at zi with strength ki and one at Z2 = z^ with strength ^2 + ^3. 

On the positive (negative) face of the plane V, or the Aq, where the normal vector to the 
face points away from (towards) the origin, we assign 7 = 1 (7 = — 1). With the same trilinear 
coordinates {xi,X2,X3), denoted by x, the point on the positive (negative) face of Ag implies 
that the vertices of Ai, or the three vortex centers zi,Z2 and Z3 in the z-plane, are in the 
counterclockwise (clockwise) orientation. Fig. [1] (a) shows the vortex centers with 7 = +1. We 
call the point x on the negative side of V the image of x on the positive side and vice versa. 
The orientations of the vortex centers or Ai can change only in crossing over an edge of Ag, 
where Ai is collinear. 

In trilinear coordinates, the two invariants ([9]) and (jlOp become, 

[k^ Xi + k2 X2 + k^ x^]p = a and (16) 

xl/''xl/''xl/'' = bp''/^'^''''^\ (17) 

With the elimination of the perimeter p, these two equations yield an integral curve, or a 
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trajectory, in V, valid for all K. It is 



'k^ k-2 h' L 1 i J J 



For a parabolic case, K = 0, Eq. ([18]) reduces to 



x^^x^x'^'^/'^^ = / or (^)'=2(^)fci=j (19) 



which is equivalent to Eq. ()17p for the parabolic case when the constant I is identified as 
Ijkik2_ jj^ ^]-^g following subsection we recount the critical points in the planar trajectories first 
presented by Synge [2] with special attention to the parabolic cases. 

2.3 Critical points R in space 

It was shown in ^ that at a critical point R, the configuration Ap is stationary, i. e., Rj = 
0, J = 1, 2, 3, if and only if Ap is either equilateral with Ri = R2 = R3 = p/3 or has zero area, 
1^1 = 0. The latter requires that the configuration Ai is coUinear lying on an edge of Aq or 
at the vertices of Aq, where two vortices coincide to one, see Eqs. dH). These critical points 
of R in space are also critical points q in plane V, or x in the trilinear coordinates, xi,X2,X3 
with p = 1, see Eqs. (|13p . An equilateral Ap, implying xi = X2 = X3 = 1/3 corresponds to 
the centroid E of Aq for positive orientation , 7 = 1 or its image E* on the opposite side for 
7 = — 1. For an equilateral configuration, it was shown that E and E* , are stable for the elliptic 
case and unstable for the hyperbolic case [2]. For the parabolic case, E and E* are degenerate 
singular points because they lye on the critical curve, which will be discussed in Sec. [3l 

For a coUinear configuration, the vortices Zj's have to lie on a straight line, C Hence we 
could use real numbers Sj to locate Zj along C, by setting si = and assigning its direction 
such that S2 = S2 — si = R3 > 0. The condition of collinearity is 

Z3- zi = v{z2 - zi), or si = 0, S2 = -R3 > 0, S3 = vs2. (20) 
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where i/ is a real number. The last equation implies R2 = ji^l-Rs- From the values of u, whether 
negative, in [0,1] or greater than 1, we see the relative positions of the vortices along £, the 
collinear configurations and the locations of the critical points on the edges of Aq (see Fig. [2] 
and [5j). For Eq. (|20p to hold for all t, we obtain from the equations of motion ([5]), that v has 
to be the root of the cubic equation, f{v) = 0, where 

f{v) = {ki + k2y - {ki + 2k2y - {ki + 2k^)v + {ki + A;3) (21) 

It was shown in [5], that there is a value K^ < 0, such that the cubic equation (j2ip has 
three unequal real roots for K > K^ and only one real root for K < K^ < in hyperbolic 
cases. For the parabolic cases, we have explicit formulas for the three roots of f{v) = 0. With 
k^ = —kik2/{ki + k2)^ we see that ki/{ki + k2) is a root of Eq. (|2ip and hence the three roots 
are: 

u = ki/{ki + k2) £ [1/2,1), Qg on side Q1Q2, (22) 



iy=[k2- Jkl + kik2 + k^]/{ki+k2) <0, Q4 on side Q2Q3, (23) 



i^ = [k2 + ^Jkl + kik2 + k^]/{ki+k2)>l, Q5 on side Q3Q1. (24) 

These three critical points lie on the respective sides of Aq. 

The trajectories defined by Eq. (jlSp and the critical points were shown in Figs. 3, 4, and 5 in 
[5] for typical elliptic , hyperbolic and parabolic cases, respectively. In this paper we study the 
critical solutions, which occur only in the parabolic case; therefore, the figure for the parabolic 
case in [5] is reproduced here in Fig. [2j 

It was noted before that a stationary point in the R space corresponds to a stationary point 
q on the plane V, but the converse may not be true. A stationary point q on P could correspond 
to similar Ap in the R space with perimeter p(t) varying. Those stationary points on V yields 
the critical curve C. We shall examine in detail the geometry of the critical curve and the planar 
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Figure 2: The trajectories for the parabolic case, kik2 + k2k^ + k^ki = 0, on the positive 
side of the plane V or the equilateral Aq. That is, the configuration Ai has counterclockwise 
orientation. Here fci = 2, A;2 = 1 and k^ = —2/3. Also shown are the coordinates, q/3, with the 
vertex Q^ as the origin and Q^E as the (3 axis. 

trajectories, in particular those intersecting the critical curve, in the next section. 

3 The planar trajectories and critical curve for the parabolic 
case 

A point q on the plane V, represented by trilinear coordinates Xj, j = 1,2,3, is related to the 
coordinates Rj^ j = 1,2,3 of the point in space by Eq. (fT3]) . We have Rj = pxj + Xjp, and 
convert Eq. ([8]) for Rj to equations for Xj and p. They are: 

xi/xi = —p/p + ^Hki{xl — xl), X2/X2 = —p/p + jHk2{xl — xl) (25) 

X3/X3 = -p/p + jHksixl - xj), where H = 2\A\/{p'^xlxlxl), (26) 
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and p/p = jH[kiXi{xl - xl) + k2X2{x1 - xl) + k3X:i{xl - xf)]. (27) 



The last equation (j27p follows either directly from Eq. ([H]) or from the sum of the preceding 
first three equations multiplied by Xj,j = 1,2,3, respectively. 
We eliminate p/p from the first three equations and obtain 

^_^ = __L7^(x)3;(x) and ^ - ^ = -l7i(x)3^(x), (28) 

Xi X2 kz Xi X2, k2 

where 3^(x) = ^2^32^1 + k^kiX2 + kik2X^. (29) 



Using the identity, J2j ^i = li 'we get J2j ^j = 0) ^^^ together with the two equations in (j28|) 
we solve for xi. We then get 

±j = -fFj n y (30) 

, j^ rX3 X2. rXi Xs X2 xi 

where Fi = xi[- — , F2 = X2[-, 7-, F3 = X3[- — . (31) 

fc2 K3 fca fei fci fc2 

To arrive at the same factor 3^(x) on the right hand sides of Eqs. (j28p and (j31|) . we used 

X = 0. Thus, the two linearly independent equations in ()3ip define the planar trajectory in V 

for parabolic cases only. These two equations plus Eq. (|27p for the variation of the perimeter 

are equivalent to the three equations dS]) for the spatial trajectory in parabolic cases. 

For points on C other than the critical points Q4, E, Q^ and E* we have the common factor 

7i > 0. For all the points on C to be stationary, it is necessary that the other common factor 

3^ = 0, 

3^(x) = k2k3xl + kiksxl + /ci/c2x| = 0, (32) 

or y(x) = -y/ks = -k2xl - kixl + {ki + k2)xl = 0. (33) 

This is the equation for the critical curve C, which exists only in the parabolic case. 

Here the derivation of the equation for C follows directly from the equations for Xj in ()3ip 
and differs from the original derivation of Synge [2j. Besides being an alternative derivation, 
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the equations in (j3ip show exphcitly the dependence of Xj's on 7 and y or Y. That is the 
direction of a trajectory is reversed when it crosses over an edge of Ag to the opposite side as 7 
changes from ±1 to ^1, respectively, or when it crosses over the critical curve C, where Y = 0, 
say from above to below as Y changes from positve to negative. These results will be needed for 
our study of the trajectories in the strip containing the critical curve in Subsecs. I3.1l and l3.2[ 
On the stationary critical curve C, the variation of the perimeter p can be obtained directly 
from one of the equations in (|25|) . say the first one. The result is 

p\t) = p\0) + A^DoSot, (34) 
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V(l-2xi)(l-2a;2)(l-2x3) 



where Dq = — ; — ; — and Sq = ; — ; — ; — -75 are constants. (35) 

K1K2 2kik2k'j,[xiX2X'j,Y 

We note that when p{t) = p(0), R is stationary. This can happen only at the centroid E 
and its image E* where xi = X2 = x^, = 1/3 and Dq = and at the critical points Q4 and Q5 
on the two sides of Ag where 1^41 = and hence 5*0 = 0. Thus the critical curve C has to pass 
through E and E* and cross over the sides Q2Q3 and Q3Q1 at the critical points Q4 and Q5, 
respectively (see Eqs. ([25]) and dMD)- 

To describe the planar trajectories and the critical curve C in the plane V, i. e., to find their 
slopes and concavities, we need to replace the trilinear coordinates satisfying the constraint, 
xi + a;2 + X3 = 1, by Cartesian coordinates in V. We choose the coordinates a and /? with the 
vertex Q3 as the origin and the bisector of /.Q3, Q3EAI, as the /? axis, where M denotes the 
midpoint of the side Q1Q2 (see Fig. [2|), with 

/3 = X3 a = {x2 — xi)/v3, with inverse transformation, (36) 

X3 = /?, xi = [1 - P - aV3]/2, and X2 = [1 - P + aV3]/2. (37) 

In terms of a and /3, Eq. (I33p for C becomes, 

3q;^ - 2\/3/ia/3 - 3/3^ + 2\/3q; - 2/? + 1 = 0, where // = (/ci - A;2)/(/ci + ^2). (38) 
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It is a hyperbola. Its intercept with the /3 axis is given by the equation 3/3^ + 2/3—1 = 0, which 
has one root, /3 = X3 = 1/3 while xi = X2 = (1 — X'i)/2 = 1/3, representing the point E or its 
image E* . The second root, [i = —1 should be excluded because P = x^ cannot be negative. 
Thus the critical curve C should be the upper branch of the hyperbola (j38p inside Aq, which 
is concave upward in the a/3 plane and passes through the centroid E (or E* on the negative 
side). 

The curve C intersects the side Q3Q1 at point x where X2 = 1/2 = xi + X3, and Eqs. ([2D|) 
and (^^ yield the equations for X3 and xi G (0,0.5), 



'kf + k^ + kik2-k2 1 ki + k2-Jkf + ki + kik2 
xs = , and xi = --X3 = . (39) 



Likewise, C intersects the side Q2Q3 at the critical point Q4, defined by Eqs. (pOj) and (j23j) . 

Thus C is divided by Q4, E, Q^ and E* into four branches; Q4E and EQ^ on the positive 
side of V and their images Q5E* and E*Qi on the negative side. From Eq. ^, we have 5*0 < 
when A ^ 0, hence p has the same sign as — 7C0, i. e., the same sign as 7(xi — X2). On 
the positive side of P, 7 = 1, we have p < on the branch EQ^ where xi < X2 < 1/2 and 
p > on the branch Q4E, where X2 < xi < 1/2. On the negative side, 7 = —1, the sign of 
p on the images of the branches is reversed. Thus the points in the branches Q4E and E*Q^, 
excluding their end points, correspond to expanding similar configurations Ap while those in 
-EQ5 and E*Q4^ correspond to contracting similar configurations. It was observed before from 
Eq. (j34p that the perimeter p remains stationary at the end points Q4,E,Q^ and E* and the 
corresponding vortex configuration Ap or R remains stationary also. 

Now we study the slopes and concavities of the trajectories, along constant I curves. Prom 
those integral curves intersecting the critical curve C, we find the variation of the invariant / 
along C and then define the boundary of the strip containing C. For this purpose, it suffices 
to study only the positive side of V, where 7 = +1. We note that the slope of a point on the 
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critical curve in the V or a/3-plane, is defined by, 

d(3/da = {kiX2 - k2Xi)V3/[2{ki + k2)x3 + /ca^i + /cixs]. (40) 

The slope has the same sign as the factor {kiX2 — k2Xi). In the af3 plane, the curve C is concave 
upward with the minimum f3 or X3 at X2 = k2Xi/ki < xi. This is in agreement with Eq. (|38|) . 
which says that C is the upper branch of a hyperbola passing through E. 
For a trajectory in the parabolic case, Eq. (fT9|) . its slope is, 

m = dp/da = V3x3{kixi - k2X2)/[2{ki + k2)xiX2 + {kixi + k2X2)xz\. (41) 

The slope m has the same sign as the factor {kixi — ^2X2), so the trajectory is concave 
downward with the maximum /5 at X2 = kixi/k2 > xi where a > 0. 

As (3 increases from to 1/2, the invariant I decreases from 00 monotonically along radial 
lines from Q3, the origin, a = 0,/3 = 0. Thus the corresponding trajectories move upward as / 
decreases. 

Since the critical curve C is concave upward while the trajectories are concave downward, C 
comes in contact with the first trajectory when they are tangent to each other. This happens 
at point E on C with a trajectory with 7=1 with common slope 

df3/da = {ki - k2)/[{ki + ^2)^3] > 0. (42) 

That is, the trajectories with J > 1 remain below C and the one with / = 1 is tangent to C at 
the centroid E. This statement holds for the images of the trajectories and the critical curve 
on the negative side of V- Hence the trajectory with I = 1 on both sides of V is the lower 
boundary of a strip S containing C. On the other hand, The trajectory with 1=1 encloses all 
the trajectories with I > 1, which are the periodic orbits around the center Q3. 

As / decreases from 1, the trajectory will intersect C at two points on the left branch Q4E 
and the right branch EQ^. respectively. To find the upper boundary trajectory of the strip S, 
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we need to find whether it is the trajectory passing through Q4 or Q5 by comparing the values 
of their invariants /(Q4) and /(Qs). For this comparison, we use Ij to denote I{Qj). Let the 
X3 coordinates of Q4 and Q5 be denoted by (3^ and /^s , then we have 

h = /(Q4) = (2/?4)-('^+')(l - 2/3^f\ and h = /(Q5) = (2/35)-('^+'Hl - 2/35)'\ (43) 

with A;2 = 1 according to Eq. ([3]). From the x^'s for Q4 and Qs, Eq. (j39p . we have 

2/34 = (A;i + l)/(g + A:i) and 2/35 = (A;i + l)/((7 + 1), (44) 



where g=-i/fc^ + A;i + l> -v/S and the ratio, 

/(A;i) = h/h = [{q - l)/{q + l)]"' [{q + k)/{q - k)] > 0. (45) 

When ki = k2 = 1, the integral curves are symmetric with respect to the /3-axis, the bisector of 
ZQa- In particular, we have /^s = /?4, / = 1 or 14 = Is, and the trajectory through Q4 coincides 
with that through Q5. 

With ki > k2 = 1, we shah show the ratio, f{ki) = h/h > 1 for ki > 1. With /(I) = 1, 
we need only prove df/dki > 0, or rather dlnf/dki > 0. With 2qdq/dk = 2ki + 1, we have 

din f/dki = - ln[(l + u)/{l -u) + 2u + u{kju'^)/{l - u^), (46) 

where v = 1/q. As ki increases from 1 to co, q increases from yS to 00, u decreases from 
l/vS to and (1 — z^^)~^ decreases from 3/2 to 1. The curve y = (1 — x^)^^ is concave 
upward and intersects the straight line or the chord, y"*" = 1 + [xu/{l — i^^)], at points (0, 1) and 
{u, (1 — u^)^^). That is, the curve lies below the chord, y'^ > y, for x G (0, v). Thus we have 

2-1/2 I r ,, 1 f' dx 1 1 + z/ ,^^, 

— — = - / y+dx > - / TT = — In (47) 

2(l-z/2) zyjo i^ Jo l-x^ 2i^ 1-iy ^ ^ 

Equation (li6]l becomes 



dlnf/dki = [(2 - i^^)u/il - u^) - ln[(l + v)/{l - v)] + {kl - l)z/V(l " ^^) > 0,(48) 
and hence 1(^4) > KQb) or Zi > -^5, for ki>l (49) 
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The equality sign holds only when ki = 1 = k2- This completes the proof, which was omitted 
in the 1988 paper |5j. Note that the trajectory through Q4 lies in the strip S, above the lower 
trajectory through E and below the upper one through Q^ as shown in Fig. [2j 

For the trajectories with I £ (Jg, Is), they will remain above C and cross over to the opposite 
side of V to form periodic orbits. For I = Ig, the trajectory becomes the separatrix through 
the critical point Qq on the side Q1Q2, as shown in Fig. [21 For I £ (0, Iq) the trajectories will 
either be periodic orbits around the center Qi or the center Q2. This completes the description 
of the trajectories, separatices and the critical curve in a parabolic case shown in Fig. [2l which 
was presented in [5]. 

From Eqs. (j49p and (j42p . we see that on the positive side of V , the critical curve C is 
contained in a strip S, bounded above by the trajectory through point Q5 on edge Q3Q1 with 
I = I5 and below by the trajectory through the centroids, E and E*, with 1 = 1. Thus the strip 
S contains besides C all trajectories with / G [-^5,1]- Likewise, we define the image of the strip 
on the negative side ofV. 

In the next subsection we will find out whether a trajectory through a point near C will be 
attracted to a point q on C or depart from q with increasing time. 

3.1 Trajectories in the strip S containing C 

Note that a point on the critical curve C is a stationary point. A trajectory will either stay 
above C or below C but will not cross over C. Thus C partitions the strip S on the positive 
side of V into two strips 5"*" and S~ above and below C. To find out whether a trajectory 
through a point in the neighborhood of C will depart from or be attracted to C, we need to 
know the direction along the trajectory as t increases. For this, we use Eqs. (|3ip . We see that 
xi and X2 changes signs only when crossing over an edge where 7 changes sign or over C where 
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y does. The factor Fi > and F2 < so that they contribute to a factor in a < because it 
has the same sign as that F2 — -Fi . That is the a will have the same sign as X2 and change its 
sign in crossing over an edge or over C. For X3 or /?, we see there is an addition sign change 
because the factor F3 can cross over zero. This sign change in /? defines the concavity of the 
trajectory, concaving downward, and the maximum of /?. Thus the directivity of the trajectories 
are defined by Eq. ([5T|) . For j = 3 we obtain 

• ,^2X2 — k^X^ . , . 

/3 = X3 = inyx^ [ uu y (50) 

K1K2 

and (3 or X3 reaches its maximum when ^2X2 = kiXi as expected from Eq. (|4ip for the slope 
TTi = 0. From the concavities of the trajectories and the critical curve, we showed in the 
preceding subsection that a trajectory with / G [/4, 1) will intersect C twice in S^ on the 
positive side of V and cross over the two edges in S- at the expanding branch, Q^E to the left 
of the centroid E, where xi > X2 with slope m > Q and at the contracting branch EQ^ to the 
right of E where X2 > Xi with ?tt, < 0. A trajectory with / G [-^5, Z4) will intersect only the right 
branch EQ^ of C and cross over the left edge in S^ and the right edge in S~ . From Eq. (j3ip . 
we see that on the positive side of P, 7 = 1, and (3 has the same sign as y[kixi — k2X2] or the 
opposite sign of my since TC > 0. 

Recalling that y > when the point is above C, and y < when below, we have the 
direction of a trajectory near C on the positive side of V, 

Near the branch EQ^, m > 0, 
Near the branch Q4,E, m < 0, 



p>o 


above C, 


(51) 


/3<0, 


below C, 


(52) 


p <0 


above C, 


(53) 


/?>0, 


below C. 


(54) 



On the negative side of P, 7 = — 1, the inequality signs will be reversed. The directions of the 
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trajectories shown in Fig. [2] are in agreement with the statements above. These statements say 

that 

In the neighborhood of a contracting (expanding) branch, EQ^ or Q^E* (Q^E or E*Q^), a 

trajectory will depart from (be attracted to) C. 

It is equivalent to say that 
When the initial vortex configuration similar to Ai or point q on a contracting branch of the 
critical curve C is disturbed off C, the trajectory will depart from C and finally be attracted to 
an expanding branch of C. 

In particular, Contracting similar vortex configuration, Ap leading to coelescence of the three 
vortices is a theoretical solution but unlikely to happen. 

Knowing the directivities of the trajectories, we shall describe in the following subsection 
how the trajectories in the strip S depart from a contracting or unstable branch of C ending on 
an expanding or stable branch and classify the trajectories into three types. 

3.2 Departure from a contracting branch of the critical curve 

For the region S^ , the centroid E divides it into two with X2 > x\ and X2 < xi respectively. 
Likewise the image of 5_ is divided by E* into two. We use 5^ to denote the subregion of S~ 
where X2 > xi, bounded above by the contracting branch EQ^ and below by the trajactory 
with I = 1 on the positive side of V and is connected across the edge where X2 = 0.5 to its 
image on the negative side of V in which the upper boundary is the expanding branch Q^E*. 
A trajectory in Sf^ with / G (-^5,1) will depart from a point below the contracting branch 
EQ^ cross over to edge and be attracted to the expanding edge E*Q^. Likewise, we use 5^ 
to denote the corresponding region where 2:2 < xi. The trajectories Sf^ have / £ (I4, 1) would 
depart from the contracting branch E*Qi and be attracted to the expanding branch Qi^E. The 
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trajactories in C~ are classified as trajectories of type I. 

For the region S^ above C, Eq. (|49p says that there is a separatrix from Q4 which intersects 
C at Si and partitions the branch EQ^ into ES4 and S4Q5 and also the strip 5+ into two strips 
S^ and S^ with / G (Zi, 1) and G {li,!^), respectively (see Fig. [2]). 

For a / G (Zi, 1), the trajectory in ^4 on the positive side of V intersects C at two points 
on the left branch Q4E and the right branch ES^, respectively, i. e., it remains on the positive 
side. The same holds for its image on the negative side. The trajectory in S^ and its image 
are of type II (see Fig. [2]). 

For a / G (15,14), the trajectory in S^ in the positive part of V will be above C, cross over 
the edge Q2Q3 and continue along the image of the trajectory in the reverse direction. The 
trajectories in ^5" are classified as type III (see Fig. [3]). 



w V S u Q, Q, u* S, V* 

A \ 1 1 / / 4 




w* Q„ 




Figure 3: Departure from a point, u or u~^ , on an unstable branch EQ^ of the critical curve 
along the trajectory u — v — w in the plane V, with 7 = +1 and along its image w — v* — u* — Uc 
in V with 7 = — 1. The last point Uc denotes the point where the trajectory intersects the 
stable branch E*Q^ of C. 

The above theoretical predictions and classifications of the trajectories in the strip S will 
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be illustrated or confirmed by the numerical examples in the next section. 

4 Numerical Simulations 

In this section, we present numerical examples to simulate or substantiate the conclusions 
arrived at in Subsec. [3.21 via the directivity of the trajectories intersecting the critical curve C. 
That is, a configuration or a point x in the neighborhold of a contracting branch of of C will 
depart from C along a trajectory with constant integral invariant / and end on an expanding 
branch of C. Instead of integrating the govening equations dH) for the vortex configurations from 
the initial configuration, the numerical examples are based on the integration of the primary 
equations of motion of the three vortices, Eqs. (j3]) in the z-plane with initial data 2;j(0)'s 
corresponding to the initial configuration Ai's in the neighborhood of C. 

First we prescribe the strengths of the three vortices belonging to the parabolic case: with 
A;2 = 1 as the scale for the strength, we choose ki = 2 and get k^ = kik2/{ki + ^2) = —2/3. 
Without losing generality, we choose the initial perimeter of Ap to be the length scale, i. e., 
p(0) = 1 and hence Rj{0) = Xj{0). 

To carry out the integration of the primary equations ^ , we shall relate the initial positions 
of the vortice4j^;j(0)'s to the initial configuration, Rj{0) = Xj{0) by (i) aligning zi and Z2 parallel 
to the imaginary axis, i. e., setting ^zi = ^Z2, (ii) requiring Qzi > Qz2 and ^^2:3 > JRzi, so 
that the configuration Ai is oriented counterclockwise, 7 = +1, and (iii) putting the stationary 
center of vorticity at the origin, z = 0. The relationships for the initial data with Rj{0) = Xj{0) 



are: 



\z3-zi\^ = Rl \Z2-Z3\^ = RI (55) 

zi — Z2 = iRs and kizi + k2Z2 + k^z^ = 0. (56) 



^Note that Xj, j = 1, 2, 3 denote the triUnear coordinates and should not be confused with ^Zj 
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Case 


Ri 


R2 


Rs 


y(o) 


7(0) 


r_ 


0.18195 


0.44396 


0.37409 


-0.00498 


0.68503 


r+ 


0.19108 


0.43424 


0.37468 


0.00501 


0.68500 


u_ 


0.10442 


0.49225 


0.40333 


-0.00500 


0.38563 


u+ 


0.10839 


0.48643 


0.40518 


0.00502 


0.38555 



Table 1: Specification of numerical examples. In all cases, ki = 2, k2 = 1, and k^ = —2/3. 



From Eq. (j55|) and the imaginary parts of the equations in (|56|) . we get 9zj(0) 



Qz2 = [-ksiRJ - Rl) + Rl{2ki + k^)]/[2R^{ki + ^2 + A^g)], 
Sj^i = 9^2 + 3^3 and Qza = — [fciQzi -|- k2'^Z2]/k^ 



(57) 
(58) 



By equating the real parts of the equations in ([56|) , we get 



3f?zi = 3f?Z2, ^z^-^zi = ^Rl-i'^sz^-^ziY and (A:i + A;2)3f?zi + fes^fiza = 0. (59) 

These three equations in turn define 3f?Zj(0). With those initial data, we carry out numerical 
integration of the primary equations (j4j), obtain the motions of the vortices, Zj{t)^s, and recover 
the vortex configuration, the sides Rj{t) and the perimeter p{t) of Ap(t), and then the three 
sides, Xj{t) = Rj/p of Ai(t). 

Ten numerical examples were carried out and the results concurred with the theoretical 
conclusions that a configuration slightly perturbed from a contracting branch of the critical 
curve C will depart from C and end on an expanding branch of C belonging to one of the three 
types described in Subsec. 13.21 Here we shall present only four typical numerical examples. 
Table 1 gives the four sets of initial configurations, Rj{0) = Xj{0) with p(0) = 1 and the 
corresponding values of 1^(0) and the invariant /(O) defined by Eqs (1331) and ([TH]) . respectively. 

Since the critical curve C is defined by Y{xj) = 0, see Eq. ([29]) . the deviation of point Xj 
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from C is measured by Y{xj) with Y positive (negative) for Xj lying above (below) the critical 
curve. 

The contracting branch EQ^ of C on the positive side of the plane V is composed of segments 
ES4 and S4Q5 with invariant / decreasing respectively from I{E) = 1 to /4 and from /4 to I5. 
See Figs. Eland El 

The integral invariants 7(0) 's in the table define the locations of the critical points on the 
contracting branch of C and each shall remain constant for t > along the trajectory towards 
an expanding branch of C. 

The first and second rows of data in Table 1 list the initial data corresponding to point r^ 
perturbed below and above, respectively, from point r in the segment ii'S'4 of C. Likewise, the 
third and fourth rows list the initial data corresponding to points u^ perturbed from a point u 
in the segment 54(55 of Cc 

For each case, the integral invariant I{t) differs from its initial value 7(0) by less than 0.1%, 
that is, the numerical solution Zj{t) does yield a planar trajectory of constant /(O). 

Figures mo and [6l show, respectively, the variations of the perimeter, log p{t), the deviation 
from the critical curve, Y{t), and the vortex configuration scaled by its perimeter, Ai, along 
the trajectories for the four cases listed in Table [H 

When a trajectory crosses over an edge to the opposite side of "P, i. e., 7 changes from ±1 to 
=pl, it continues along the image of that on the positive side but in reverse direction; therefore, 
at an edge, logp reaches an extremum and so does Y while the three vortices become collinear. 
We also note that on either side of P or the a/3 plane, the trajectories, 7 = constant, are concave 
downward while the constant Y{xj) lines are concave upward similar to that of C where Y = 0. 



^The Xj^s of point r on C, where Y — can be defined from the hnear interpolation of Sj's of r with respect 
to Y (0). Likewise we can locate the Xjis of point u on C from those of u . 
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Figure 4: Variation of Y(x) along a trajectory, x(t), departing from an unstable branch of the 
critical curve to a stable branch. Point of departure: r^ , r+, u^ , and n+ in a), b), c) and d). 

The trajectory 1 = 1 and the critical curve 1" = are of opposite concavity. They are tangent 
to each other a E (or E*). The trajectories intersecting the left (right) branch Q4E or EQ^ of 
C will have I{Qa) < / < 1 or liQ^) < I < 1. This statement holds also for the images of the 
trajectories. 

The trajectories with / < 1 will be tangent to a constant y-line with Y > above C. Thus 
along a trajectory with / < 1, ^(x) reaches a maximum at a point above C. Again this holds 
also for the image trajectory on the negative side of V 

Along a trajectory initiating at point below C, say r~ or n^, the deviation, 5 = —Y, (the 
perimeter or —logp) will increase from almost zero, or 5 ^ 1, to a maximum at the edge and 
then return to zero as the trajectory ending at a point on the stable branch of C while Ai(t) 
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Figure 5: Variation of logp(x) along a trajectory, x(t), departing from an unstable branch of 
the critical curve to a stable branch. Point of departure: r^, r+, u^ , and u^ in a), b), c) and 
d). 

becomes equivalent to Ai(0) with a change of orientation. The trajectory is of type I. See part 
a) and c) of these three figures. 

For a trajectory initiating from a point, say r"*", above the unstable branch ES4 of C where 
I{Si) = I{Qi) > liQ^), the deviation 6 will increase from a small initial value , <C 1, to a 
maximum and then return to zero as the trajectory arrives at a point on the stable branch Q4^E 
of C without crossing over an edge, i. e., staying on the same side of V, as in Fig. [Hb). The 
trajectory is of type II. The perimeter will decrease from 1 to a minimum and then keep on 
increasing as t increases. See part b) in Fig. [5j Note that the final configuration corresponding 
to the point on the left branch Q4E is different from the initial one on the right branch £'^4, 
as shown in Fig. [6] d) . 
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Figure 6: Evolution of Ai departing from the unstable branch on the positive side of V to the 
stable branch on the negative side of type I or III, or to that on the positive side of type II. 

For a trajectory initiated from a point, u^ above C with I G {I{Qb)^I{QA))-, the deviation, 
5 = ^(0), will increase from a small initial value ^ 1 to a maximum, crossing the edge Q2Q3 
with Y reaching a local minumum and continue to the negative side of V with Y appearing 
as the image of the first part passing through the image point of u'^ , then return to zero as it 
reaches a point on the stable branch of C. The image of the point will be close to the initial 
point u^. The trajectory is of type III. 

Our numerical examples show the three distinct types of trajectories expected in the the- 
oretical studies in the preceding section. Note that the deviation Y{t) and the perimeter p{t) 
for type III have three extrema as t increases while those for type I and II have only one. For 
type I and III, there are changes of orientation and the final configurations are nearly similar 
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to the original ones with perimeter increasing. For type II, there are no changes of orientation 
but the final expanding configuration is not similar to the initial one. 

The trajectories depart from the unstable branch E*Q4 of C on the negative side of V. 
Because of I4 > I^, the departing trajectories are only of two types: trajectories of type I begin 
from a point below the branch Q4E* while those of type II begin from a point above. 

5 Conclusion 

This paper presents a complete description of the grobal dynamics of three point vortices 
initially in the neighborhood of the planar critical curve C in trilinear coordinates Xj, j = 
1,2,3. The trilinear coordinates, only two of which are independent because of the constraint 
xi + X2 + X3 = 1, were introduced by Synge (1949) to study the deformation of the triangle Ap 
formed by the three vortices with sides Rj and perimeter p via that of a similar triangle with 
sides Xj and perimeter one. A spatial trajectory Rj{t) is then reduced to a planar trajectory 
Xj{t). In the plane, a stationary critical curve C was found by Synge (1949) for parabolic cases. 
Each point on C corresponds to a radial spatial trajectory either moving away or towards the 
origin with similar Ap expanding p > 0, or contracting p < 0, respectively. The latter would 
lead to coelescence of three vortices. According to the spatial trajectories, the critical curve 
is partitioned into expanding or contracting branches. It was shown by Tavantzis and Ting 
(1988) that the expanding radial (spatial) trajectories are stable while the contracting ones 
are unstable. In this paper we show that a contracting radial trajectory slightly disturbed will 
depart from the radial line, move along a spatial trajectory nonsimilar to the initial configuration 
and eventually approach a radial expanding trajectory. Thus we come to the conclusion that 
three distinct vortices almost never coelescence. 
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